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Summary 

Frangible  Armour  Piercing  Discarding  Sabot  (FAPDS)  munition  is  a  new  type  of  medium  calibre 
munition  for  air  defense  systems.  Its  penetration  behaviour  has  been  studied  by  TNO-PML  in 
cooperation  with  BWB,  IABG,  and  Rheinmetali.  These  studies  lead  to  a  simulation  model 
describing  the  penetration  behaviour  of  35  mm  FAPDS  projectiles  in  thin  plate  arrays. 

This  report  discusses  an  improved  version  of  the  semi-empirical  terminal  ballistics  model  for 
FAPDS  munition.  With  regard  to  the  first  simulation  model  (as  reported  in  report  PML  1992-6), 
the  characteristics  of  the  fragmentation  process  have  been  changed  (use  of  Mott  distribution 
instead  of  random  distribution,  increase  in  number  of  particles,  and  decrease  in  scattering  angle), 
and  the  parasitic  effect  has  been  described  more  accurately  (effect  of  hole  and  particle  size  on 
penetration  behaviour  is  included).  Due  to  the  changes  in  the  model,  the  simulation  results  are 
consistent  with  the  experimental  results.  Now,  the  maximum  penetration  capacity  and  the 
fragment  beam  width  can  be  calculated  accurately  for  35  mm  FAPDS  munition.  The  influence  of 
the  target  array  composition  and  the  rifling  twist  of  the  gun  on  the  simulation  results  can  also  be 
predicted  well. 

However,  some  problems  still  have  be  solved  before  the  terminal  ballistics  model  can  be 
incorporated  into  a  vulnerability  assessment  code.  The  major  prooiem  is  finding  a  way  to  deal  with 
synergistic  effects.  It  is  expected  that  the  problems  will  be  solved  within  the  next  half  year,  so  that 
effectiveness  calculations  can  be  performed  against  aircraft  and  helicopters.  The  calculation  results 
will  be  used  for  the  evaluation  of  several  munition  types  as  replacement  of  todays  munition  for  air 
defense  systems. 

Samenvatting 

Frangible  Armour  Piercing  Discarding  Sabot  (FAPDS)  munitie  is  een  nieuw  type  middenkaliber 
munitie  voor  luchtverdedigingssystemen.  Het  penetratiegedrag  van  dit  type  munitie  is  door  PML- 
TNO  in  samenwerking  met  BWB,  IABG  en  Rheinmetali  bestudeerd.  Deze  studies  hebben 
geresulteerd  in  een  simulatiemodel  dat  het  penetratiegedrag  van  35  mm  FAPDS  projectielen  in 
dunne  platenarrays  beschrijft. 

Dit  rapport  geeft  een  beschrijving  van  een  verbeterde  versie  van  het  semi-empirische 
eindballistische  model  voor  FAPDS  munitie.  Ten  opzichte  van  de  eerste  versie  van  het 
simulatiemodel  (zie  rapport  PML  1992-6)  zijn  in  de  verbeterde  versie  de  eigenschappen  van  het 
verschervingsproces  aangepast  (gebruik  van  een  Mott  verdeling  in  plaats  van  een  random 


I 


verdeling,  grotere  aantallen  scherven  en  een  kleinere  spreidingshoek)  en  is  de  beschrijving  van  het 
parasitaire  effect  verbeterd  (het  effect  van  de  scherf-  en  gatgrootte  op  het  penetratiegedrag  wordt 
in  rekening  gebracht).  Door  deze  veranderingen  komen  de  simuiatieresuitaten  overeen  met  de 
experimentele  resultaten.  De  maximale  penetratiecapaciteit  en  de  bundelbreedte  kunnen  voor  35 
mm  FAPDS  munide  nauwkeurig  worden  berekend.  Bovendien  zijn  de  invloed  van  de 
doelsamenstelling  en  van  de  spoed  van  het  kanon  op  de  resultaten  goed  te  beschrijven  met  het 
model. 

Er  zijn  echter  nog  enkele  problemen  die  de  implemenratie  in  een  kwetsbaarheidscode 
bemoeilijken.  Het  grootste  probleem  wordt  veroorzaakt  door  het  synergetisch  effect.  De 
problemen  zullen  naar  verwachting  binnen  een  half  jaar  worden  opgelost,  zodat 
effectiviteitsberekeningen  tegen  vliegtuigen  en  helicopters  kunnen  worden  uitgevoerd.  De 
resultaten  van  deze  berekeningen  zullen  worden  gebruikt  bij  de  evaluatie  van  verschillende  typen 
munide  voor  de  luchtverdediging. 
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1  INTRODUCTION 

One  of  the  new  medium  calibre  munition  types  is  the  so-called  Frangible  Armour  Piercing 
Discarding  Sabot  (FAPDS)  munition.  The  FAPDS  round  consists  of  a  sub-calibre  projectile, 
which  may  comprise  a  frangible  core  enclosed  by  a  steel  nose  and  a  tungsten  rear  part.  The 
frangible  part  is  meant  to  continuously  break  up  during  penetration  of  the  target  and  to  generate  a 
fragment  beam.  The  fragment  beam  gives  an  axial  as  well  as  a  radial  damage  pattern. 

FAPDS  munition  has,  as  a  sub-calibre  projectile,  a  low  time-of-flight  and  therefore  relatively  high 
hit  probability  and,  because  of  its  design,  can  defeat  different  kinds  of  targets,  i.e.  soft  material 
targets  such  as  aircraft,  and  lightly  armoured  targets  such  as  helicopters  and  personnel  carriers. 

The  Bundesamt  fur  Wehrtechnik  und  Beschaffung  (BWB)  and  the  Rheinmetall  company  have 
been  performing  experiments  to  study  the  penetration  behaviour  of  the  FAPDS  projectile. 
However,  no  satisfactory  description  has  as  yet  been  found  to  describe  the  penetration  process. 
Such  a  physical  description  is  necessary  to  develop  a  penetration  model  that  can  be  incorporated 
into  vulnerability  assessment  codes. 

In  1986,  TNO-PML  was  asked  by  Directie  Materieel  Koninklijke  Landmacht,  afdeling  Munitie 
(DMKL)  to  contribute  to  the  development  of  a  physical  description  of  the  FAPDS  projectile 
behaviour.  This  was  done  in  cooperation  with  the  BWB,  the  Industrieanlagen  Betriebsgesellschaft 
(IABG)  and  Rheinmetall. 

Since  1987,  TNO  has  been  performing  an  extensive  experimental  programme  [1,2]  to  study  the 
terminal  ballistics  of  the  FAPDS  munition.  The  study  has  concentrated  on  a  concept  of 
Rheinmetall.  Experiments  on  multiple  plate  arrays  have  led  to  a  better  qualitative  understanding  of 
the  35  mm  FAPDS  munition  penetration  behaviour.  As  a  result,  quantifications  of  effects  involved 
in  the  penetration  process  could  be  made  [3,4]. 

In  the  second  phase  of  the  study,  started  in  June  1991,  special  attention  was  paid  to  the  modelling 
of  the  terminal  ballistics  behaviour  of  FAPDS  munition  for  vulnerability  assessment  studies.  The 
incorporation  of  such  a  model  into  vulnerability  assessment  codes  is  rather  complicated,  as 
indicated  below. 

In  ordinary  vulnerability  assessment  codes,  such  as  the  TNO-PML  TARget  Vulnerability 
Assessment  Code  (TARVAC),  the  projectile  trajectory  is  simulated  by  a  single  (straight)  line 
which  is  traced  through  the  target.  Along  this  line,  calculations  are  performed  to  determine  the 
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extent  of  the  projectile’s  penetration  into  the  target  and  the  damage  inflicted  upon  the  target.  This 
approach  leads  to  satisfactory  results  for  single  projectile  munition  types. 

In  contrast,  the  FAPDS  munition  desintegrates  into  many  fragments  or  particles1  which  penetrate 
the  target  more  or  less  simultaneously,  while  the  fragments  continuously  break  up  into  little 
particles  whenever  they  hit  some  part  of  the  target.  Therefore,  simulation  of  the  projectile 
trajectory  with  a  single  shot  line  is  not  sufficient.  A  shot  line  should  be  evaluated  for  each 
fragment.  Furthermore,  the  effects  produced  by  the  individual  fragments  are  not  independent  for 
two  reasons. 

1  Whenever  a  fragment  creates  a  hole  in  the  target,  another  fragment  which  arrives  at  a  later  time 
can  fly  straight  through  the  hole.  The  latter  does  not  need  to  perforate  the  target  again,  and 
hence  will  not  break  up.  This  effect  is  referred  to  as  parasitism. 

2  The  impact  of  multiple  fragments  can  give  a  greater  effect  than  an  equal  number  of  fragments 
hitting  the  target  separately.  In  other  words,  the  total  effect  will  be  greater  than  the  sum  of  the 
parts  due  to  synergism. 

A  first  approach  for  an  FAPDS  simulation  model  to  describe  the  penetration  behaviour  of  35  mm 
FAPDS  projectiles  into  thin  plate  arrays  was  presented  in  [5J.  In  the  mean  time,  the  FAPDS 
simulation  model  has  been  improved,  and  some  aspects  have  been  studied  in  greater  detail.  This 
report  discusses  the  improved  simulation  model. 

The  simulation  model  is  valid  for  the  FAPDS  round  under  development,  which  is  presented  in 
Chapter  2.  In  Chapter  3,  a  full  description  of  the  improved  FAPDS  simulation  model  is  presented. 
The  comparison  of  the  results  achieved  with  the  simulation  model  with  the  experimental  results  is 
discussed  in  Chapter  4,  as  well  as  the  influence  of  certain  model  parameters  on  the  simulation 
results.  Finally,  the  problems  related  to  the  implementation  of  the  simulation  model  in  a 
vulnerability  assessment  code  are  discussed  in  Chapter  5. 

It  is  expected  that  these  problems  will  be  solved  within  the  next  half  year,  so  eventually 
calculations  can  be  performed  against  aircraft  and  helicopters.  The  calculation  results  will  be  used 
for  the  evaluation  of  several  munition  types  as  replacement  of  todays  munition  for  air  defence 
systems. 
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We  will  use  the  words  "fragments"  and  "particles"  to  indicate  the  difference 
between  respectively  the  fragments  before  and  after  breaking  up. 
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2  THE  35  MM  FAPDS  PROJECTILE 


The  FAPDS  simulation  model  has  been  validated  for  the  35  mm  FAPDS  round  under 
development,  which  comprises  a  frangible  core  enclosed  by  a  steel  nose  and  a  tungsten  rear  pan. 
This  muntion  was  produced  by  Rheinmetall.  A  cross-section  of  the  projectile  is  shown  in  Figure  1 . 


Tungsten  alloy 


Breakup  material 


Figure  1  35  mm  FAPDS  projectile  (penetrator  and  sabot) 


At  this  moment,  also  different  kinds  of  FAPDS  projectiles  (e.g.  projectiles  completely  made  of 
frangible  material)  are  available.  Even  different  calibre  FAPDS  rounds  (e.g.  25  mm)  are  being 
produced.  The  FAPDS  simulation  model  has  not  been  validated  for  these  projectiles  yet. 
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3  THE  FAPDS  SIMULATION  MODEL 

On  the  basis  of  the  experimental  results  [1 ,2,3,6],  an  FAPDS  simulation  model  has  been 
developed.  This  model  describes  the  penetration  behaviour  of  35  mm  FAPDS  projectiles  into  thin 
plate  arrays. 

The  FAPDS  simulation  model  is  a  cascade  model;  the  frangible  part  of  the  projectile  will  break  up 
on  target  impact,  and  each  fragment  breaks  up  into  smaller  particles  at  each  subsequent  target 
plate  encounter.  The  perforation  process  of  a  target  plate  is  divided  into  two  successive  steps;  the 
penetration  of  the  target  plate,  and  the  breaking  up  of  a  fragment  after  perforation.  A  penetration 
model  is  used  to  determine  the  residual  mass  and  velocity  of  the  fragments  after  perforation  of  the 
target  plate.  A  fragmentation  model  is  used  to  calculate  mass,  velocity  and  spatial  distributions  for 
the  particles  into  which  a  fragment  breaks  up.  Finally,  a  model  to  account  for  the  effect  of  parasitic 
fragments  is  included  in  the  FAPDS  simulation  model. 

The  synergistic  effect  is  not  included  in  the  simulation  model.  However,  an  estimate  of  the 
synergistic  effect  will  be  made  while  evaluating  the  simulation  results  (Chapter  4). 

The  respective  models,  which  form  together  the  FAPDS  simulation  model,  are  described  below. 
The  models  have  been  implemented  in  an  FAPDS  simulation  computer  program.  This  computer 
program  is  discussed  at  the  end  of  this  chapter. 

3. 1  The  penetration  model 

In  the  FAPDS  simulation  model  the  so-called  THOR  equations  are  used  to  calculate  the  residual 
mass  and  velocity  of  the  fragments  after  penetration  of  a  target  plate.  The  equations  are  also  used 
to  describe  the  penetration  of  the  nose  and  rear  part  of  the  FAPDS  projectile.  However,  these 
parts  will  not  break  up  after  perforation  of  the  target. 

The  THOR  equations  were  derived  experimentally  from  ballistic  tests  [7]  and  will  be  discussed 
briefly. 

Consider  the  configuration  shown  in  Figure  2. 
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Figure  2  A  fragment  penetrating  a  target 

A  fragment  with  mass  mjn  and  velocity  Vjn  hits  the  target.  The  angle  between  the  velocity  vector  of 
the  fragment  and  the  normal  of  the  target  surface  is  defined  as  0.  During  the  penetration,  the 
velocity  of  the  fragment  as  well  as  its  mass  reduces.  If  the  fragment  strikes  the  target  with  sufficient 
energy,  it  will  perforate1  the  target.  The  residual  mass  mout  and  velocity  vout  can  be  calculated 
with  aid  of  the  following  THOR  equations: 

vout  =  vin  *  loCl  *  (t*A)C2  *  minC3  *|^r^jC4  *  vinC5  (31) 

mout  =  min  -  1 0K 1  *  (t*A)K2  *  minK3  *  j  K4  *  vinK5  (3.2) 

where  C1..C5  and  K1..K5  are  material  constants 
t  =  target  thickness 

A  is  the  average  presented  area  of  the  fragment  calculated  from  the  shape  factor  S 
and  the  mass  mjn: 

A  =  S  *  m;^3  (3.3) 

1  We  use  the  word  "perforate"  to  indicate  that  a  fragment  is  capable  of 

penetrating  a  target  plate  completely,  and  the  residual  mass  and  velocity  of 
the  fragment  are  greater  than  zero. 
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If,  as  result  of  equations  3. 1  and  3.2,  the  residual  mass  or  velocity  is  less  than  or  equal  to  zero,  the 
fragment  is  not  capable  of  perforating  the  target. 

Since  within  the  beam  a  velocity  gradient  is  present,  a  finite  probability  exists  that  a  fragment 
encounters  a  hole  in  a  target  plate,  created  by  a  fragment  arrived  earlier.  Depending  on  the 
dimensions  of  the  fragment  and  of  the  hole,  the  fragment  can  fly  through  that  hole  without  any 
interaction  with  the  target  plate  (parasitism). 

If  the  fragment  is  larger  than  the  hole  encountered,  it  will  interact  with  the  target  plate  and  it  will 
break  up.  We  have  derived  modified  THOR  equations  to  calculate  the  residual  mass  and  velocity 
of  the  fragments  for  this  specific  situation.  An  extensive  description,  as  well  as  the  derivation,  of 
the  modified  THOR  equations  can  be  found  in  Annex  A. 

The  principle  of  these  equations  will  be  explained  briefly.  The  equations  are  based  on  energy 
conservation  laws.  First  of  all,  we  calculate  the  minimum  energy  necessary  to  create  the  hole  in  the 
target  plate.  This  minimum  energy  is  added  to  the  residual  energy  of  the  fragment  necessary  to 
perforate  the  same  target  plate  without  a  hole.  From  the  resulting  energy,  the  residual  mass  and 
velocity  of  the  fragment  is  deduced.  The  original  THOR  equations  (3.1  and  3.2)  are  applied  to 
calculate  the  minimum  energy,  and  to  deduce  the  residual  mass  and  velocity. 

An  important  restriction  of  the  modified  THOR  equations  is  that  they  can  be  used  only  for  small 
impact  angles.  However,  this  is  hardly  a  problem  for  FAPDS  simulations,  because  the  fragments 
are  limited  to  a  narrow  cone  (Chapter  4). 

Another  important  parameter  needed  for  the  FAPDS  simulation  is  the  perforation  time  At.  The 
perforation  time  represents  the  time  it  will  take  a  fragment  to  perforate  the  target.  The  principle  of 
constant  force  is  applied  in  the  simulation  model  to  calculate  the  perforation  time.  On  the  basis  of 
the  momentum  and  energy  conservation  laws  (3.4  and  3.5),  we  arrive  at  the  following  expression 
(3.6)  for  the  perforation  time: 


At 

min  *  vin  ~  mout  *  vout  =  =  ^  *  At 


(3.4) 


1/2  mjn  *  vjn^  -  1/2  mout  *  vout^  - 


L 

fF*ds  =  F  *  L 


(3.5) 
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L  *  (m,n  *  vjn  -  mout  *  vout) 

At  = - — - - >i - t  (3.6) 

(1/2  mjn  *  Vjn  -  1/2  mout  *  vout  ) 

where  F  =  Force 

L  =  Line-of-sight  thickness 

3.2  The  fragmentation  model 

Because  quantitative  theories  governing  the  fragmentation  process  have  not  been  developed  yet, 
we  have  used  experimental  results  to  obtain  the  spatial,  the  fragment  mass  and  the  velocity 
distributions.  In  addition,  experiments  with  frangible  Fragment  Simulating  Projectiles  (FSPs)  are 
planned  to  obtain  more  information  on  the  fragmentation  of  individual  fragments. 

The  elements  of  the  fragmentation  model  are  described  below. 

Mass  distribution 

We  have  assumed  that  a  single  fragment  (mass  m)  breaks  up  according  to  a  Mott  distribution  into 
Nq  smaller  particles  (masses  m,,  i  6  [1,Nq]).  The  Mott  distribution  is  a  cumulative  distribution, 
and  is  defined  by  equation  3.7. 


N(mf )  =  N0  *  exp-  (  |  2  (3-7) 

v  0.5*  mV 

where  N(rrij)  =  cumulative  number  of  fragments  with  a  mass  larger  than  or  equal  to  mj 
m  =  average  fragment  mass  (=  m  /  Nq  ) 

The  total  mass  of  the  particles  should  be  equal  to  the  original  fragment  mass: 


Experiments  [2]  indicate  that  the  number  of  particles  into  which  a  fragment  breaks  up,  ranges 
from  10  to  20. 
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Velocity  and  spatial  distribution 

After  breaking  up,  the  resulting  particles  travel  with  different  velocities  ( vj,  i  e  [1,Nq])>  both  in 
magnitude  and  in  direction.  In  the  simulation  model,  the  magnitude  of  the  velocity  has  been 
chosen  equal  for  all  particles  originating  from  one  fragment.  The  velocity  of  the  particles  is  equal 
to  the  residual  velocity  of  the  fragment  as  calculated  with  the  (modified)  THOR  equations. 

The  difference  in  direction  is  accounted  for  by  a  scattering  angle  aji 


COSCtj 


Vj*V 

fv;|*|v| 


(3.9) 


where  v  is  the  velocity  vector  of  the  original  fragment  before  break  up. 

The  scattering  of  the  individual  particles  is  limited  to  a  narrow  cone  (aperture  angle  about 
3  degrees)  [2],  The  simulation  model  uses  a  normal  distribution  for  the  scattering  of  the  panicles 
within  this  cone. 


The  rifling  twist  of  the  gun,  from  which  the  FAPDS  round  is  fired,  causes  the  projectile  to  spin. 
The  spin  rate  has  great  influence  on  the  spatial  scattering  of  the  particles  [6] .  The  spin  rate  (to)  of 
the  projectile  (diameter  D)  depends  on  the  rifling  twist  (RT),  the  calibre  C,  and  the  muzzle 
velocity  vg  according  to: 


tan(RT) 

““1C72T 


* 


v0 


(3.10) 


A  particle  which  breaks  off  from  the  frangible  part  of  the  original  projectile  can  be  imparted  a 
maximum  tangential  velocity  (vt  max)  equal  to: 


vt  max 


=  tan(RT)  *  vq  *  ^ 


(3.11) 
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At  random,  a  fraction  of  the  maximum  tangential  velocity  is  added  to  the  particle  velocity  in  the 
simulation  model.  For  particles  originating  from  the  original  projectile,  a  first  order  estimate  can 
be  made  for  the  contribution  of  the  tangential  velocity  to  the  scattering.  The  maximum  scattering 
angle  (amax)  is  given  by  equation  (3.12). 


,  ,  vt.max  *  D 

tan(«max)  -  =  tan(RT)  ^ 


(3.12) 


amax  ~  RT  *  ^  ,  for  small  values  of  RT  (3.13) 

Equation  (3.13)  shows  that  an  increase  in  rifling  twist  will  lead  to  a  proportional  increase  in 
scattering  angle.  As  a  result,  the  fragment  beam  will  become  wider. 

For  the  35  mm  FAPDS  round  under  study,  the  following  data  have  been  used  [6] : 

D  =18  mm 

Calibre  =  35  mm 
RT  =  6' 

As  a  result,  the  maximum  scattering  angle  (amax)  will  be  about  three  degrees. 

3.3  The  parasitism  model 

As  stated  before,  penetrating  fragments  create  holes  which  will  allow  other  fragments  to  fly 
through  the  target  without  interactions.  The  simulation  model  detects  whether  a  fragment  flies 
through  an  existing  hole.  If  the  fragment  is  smaller  than  the  hole,  it  will  fly  through  that  hole.  If  the 
fragment  is  larger  than  the  hole,  the  fragment  is  not  classed  as  being  parasitic.  Hence,  it  will  break 
up  after  perforation. 

3.4  The  FAPDS  simulation  computer  program 

The  models  discussed  have  been  implemented  in  an  FAPDS  simulation  computer  program.  The 
program  is  capable  of  performing  a  complete  simulation  loop,  given  a  normal  impact  of  a  35  mm 
FAPDS  projectile  on  a  thin  plate  array.  As  a  result  of  the  simulation,  a  list  with  holes  is  generated 
which  yields  the  damage  inflicted  upon  the  target.  Furthermore,  statistical  information  on  the 


TNO-report 


PML  232491  308 


Page 

14 


fragments  can  be  evaluated  (e.g.  mass  distribution).  The  information  will  be  used  to  compare  the 
simulation  results  with  the  experimental  results  (Chapter  4). 

Two  major  mechanisms  determine  the  structure  of  the  FAPDS  simulation  computer  program:  the 
breaking  up  of  the  fragments  after  each  perforation  and  the  effect  of  parasitic  fragments. 

The  breaking  up  is  simulated  by  creating  a  new  set  of  vectors  after  each  penetration.  Each  vector 
represents  a  single  shot  line.  The  intersections  of  the  shot  lines  with  the  target  (geometrical  model) 
are  calculated  with  the  aid  of  ray-tracing  techniques. 

The  effect  of  parasitic  fragments  implicates  a  time  relation  among  the  fragments.  As  a  result,  the 
time-element  must  be  incorporated  in  the  simulation.  We  have  chosen  to  make  the  simulation 
event-driven.  In  this  approach,  each  impact  of  a  fragment  on  the  target  is  regarded  as  an  event.  All 
events  are  stored  in  an  event  list.  The  first  event  in  the  simulation  loop  is  the  impact  of  the  original 
projectile.  The  event  list  is  empty  if  no  fragments,  which  are  capable  of  perforating  the  target,  are 
left.  Then,  all  fragments  are  stopped  within  the  target  or  have  left  the  plate  array  at  a  side. 

The  parasitism  model  itself  is  implemented  through  the  maintenance  of  a  list  with  holes.  A  new 
hole  is  added  to  the  list  each  time  a  fragment  perforates  the  target.  A  hole  is  simply  described  by  a 
hollow  cylinder  in  the  target.  The  radius  r  of  the  cylinder  is  calculated  from  the  average  presented 
area  (A)  of  the  fragment: 


A  =  S  *  m2/3 


(3.14) 


and 


A  =  n  *  r2 

Combining  equations  3.14  and  3.15,  the  radius  of  the  cylinder  equals: 


(3.15) 


(3.16) 


To  determine  whether  a  fragment  flies  through  a  hole,  the  list  with  holes  is  scanned  until  a  hole  is 
found  for  which  the  impact  location  on  the  target  is  located  within  the  cylinder.  If  the  presented 
area  of  the  fragment  is  larger  than  the  hole  area,  the  hole  size  will  be  adapted  if  the  modified 
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THOR  equations  predict  that  the  relevant  fragment  is  capable  of  perforating  the  target.  If  no  hole 
is  found,  a  new  hole  is  added  to  the  list  for  that  fragment  if  the  relevant  fragment  is  capable  of 
perforating  the  target. 

The  complete  structure  of  the  FAPDS  simulation  program  and  some  details  on  the  evaluation  of 
individual  events  are  described  in  Annex  B. 

The  main  assumptions  and  restrictions  of  the  computer  program  regarding  the  various  models  are 
summarized  below. 

Penetration  model: 

-  The  (modified)  THOR  equations  can  be  used  to  describe  the  penetration  process  for  all 
fragments,  and  for  the  main  parts  of  the  FAPDS  projectile. 

-  The  FAPDS  projectile  strikes  the  target  (thin  plate  array)  perpendicularly. 

-  The  shape  factor  is  equal  for  all  fragments. 

-  The  effect  of  tumbling  is  ignored,  as  well  as  air  drag. 

-  The  hole  size  is  equal  to  the  presented  area  of  the  fragment  that  created  the  hole. 

Fragmentation  model: 

-  Except  for  the  parasitic  fragments,  each  fragment  which  perforates  the  target  breaks  up  into  a 
number  of  particles  according  to  a  Mott  distribution.  The  actual  number  is  drawn  randomly 
out  of  a  range,  which  can  be  defined  by  the  user. 

-  After  a  fragment  is  broken  up,  the  magnitude  of  the  velocity  is  equal  for  all  particles. 

-  The  spatial  distribution  of  the  scattering  angle  is  normal.  The  maximum  scattering  angle  can 
be  defined  by  the  user. 

-  The  spin  of  the  original  FAPDS  projectile  is  lost  after  the  first  fragmentation. 

Overall: 

-  All  events  are  regarded  to  be  independent  (except  for  parasitism  effects),  meaning  no  collisions 
or  synergistic  effects  are  included. 
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4  VALIDATION  OF  THE  FAPDS  SIMULATION  MODEL 


To  validate  the  FAPDS  simulation  model  and  its  implementation  ,  the  results  of  the  experimental 
programme  performed  by  TNO-PML  have  been  used  [1,2].  The  experiments  were  conducted 
with  35  mm  FAPDS  projectiles  against  the  so-called  TNO-arrays.  The  TNO-array  consists  of 
aluminium  plates;  the  first  plate  is  8  mm  thick,  all  other  plates  are  2  mm  thick.  All  plates  are 
placed  normal  to  the  line-of-fire  at  a  mutual  distance  of  30  cm.  In  the  experiments  up  to  seven 
plates  were  used  for  the  TNO-array. 

The  simulation  and  the  experimental  results  are  compared  in  the  first  part  of  this  chapter.  In  the 
second  part,  we  demonstrate  the  influence  of  target  composition  and  rifling  twist  on  the  results. 
Finally,  the  influence  of  certain  parameters  of  the  simulation  model  on  the  results  is  discussed. 

4. 1  FAPDS  simulation  results 


The  FAPDS  simulation  computer  program  has  been  used  to  perform  simulations  of  35  mm 
FAPDS  projectile's  penetration  into  the  TNO-array.  The  values  of  the  relevant  parameters,  that 
have  been  used  in  the  simulation  runs,  are  listed  below. 

35  mm  FAPDS  projectile: 
mass  front  part 
mass  frangible  part 
mass  rear  part 

shape  factor  of  the  fragments/particles 
impact  velocity 
impact  angle 
rifling  twist  of  the  gun 

number  of  particles  into  which  a  fragment  breaks  up  (Nq) 
scattering  angle  for  the  individual  fragments  (aj) 

The  values  of  the  last  two  parameters  are  not  identical  to  the  experimentally  obtained  values,  but 
they  proved  to  give  better  overall  results.  Explanations  to  justify  these  values  are  presented  in  this 
chapter. 


1 1-79  g 
244.91  g 
43.59  g 

0.004288  m2  /  kg2/3 
1400  m/s 
ONATO 
6’ 

10<N0<  100 
2 • 


Target  damage 

Figure  3  shows  the  damage  caused  by  the  fragments  to  the  plates  of  the  TNO-array  as  calculated 
with  the  FAPDS  simulation  program.  This  figure  clearly  demonstrates  that  synergism  is  not 
included  in  the  simulation  model,  otherwise  several  big  holes  surrounded  by  zones  of  indents  and 
cracks  would  be  visible  in  the  target  plates. 
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Figure  3  The  damage  caused  to  the  plates  of  the  TNO-array  by  a  35  mm  FAPDS  projectile 
(the  holes  and  the  non-perforating  impacts  are  shown) 

From  the  experiments  it  was  deduced  that  the  perforation  pattern  in  a  plate  can  be  typically 
described  by  three  damage  areas  [1];  a  perforated  zone,  a  zone  which  surrounds  the  perforated  one 
and  the  non-perforating  impacts,  and  a  zone  which  encompasses  all  outward  running  cracks 
(Figure  4A).  The  hole  diameter  is  defined  as  a  characteristic  dimension  for  the  perforated  ..one, 
and  it  is  determined  by  idealizing  the  perforated  zone  to  a  circle  with  an  area  equal  to  the 
measured  hole  area. 
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Figure  4A  Definition  of  characteristic  Figure  4B  Definition  of  characteristic 

diameters  (experiments)  diameters  (according  to 

simulations) 

In  Figure  5,  the  hole  diameters  as  measured  in  the  TNO-array  and  as  calculated  by  the  FAPDS 
simulation  program  are  shown.  The  hole  diameter  according  to  the  simulations  is  defined  as  the 
diameter  of  the  smallest  circle  which  encompasses  all  holes  in  the  target  plate  (Figure  4B).  As  a 
result,  the  hole  diameters  according  to  the  simulations  will  be  larger  than  the  measured  ones.  By 
defining  the  hole  diameter  in  such  a  way,  the  synergistic  effect  has  been  introduced  in  the 
evaluation  of  the  simulation  results.  The  hole  diameter  is  an  estimate  of  the  hole  area  that  would 
be  measured  in  real  targets. 


TNO-report 


PM1. 232491308 


Page 

19 


Figure  5  Experimental  and  theoretical  hole  diameters  of  the  TNO-array 

(★  experiments,  ▲  experiments  in  which  the  hole  diameters  were  similarly  defined 
as  in  the  simulations,  X  simulations) 

The  definition  of  the  hole  diameter  as  applied  in  the  simulation  results  has  also  been  applied  to  the 
experimental  results.  An  estimate  of  this  hole  diameter  was  made  bv  measuring  the  impact  pattern 
in  the  fragment  catcher,  which  was  placed  behind  the  last  plate  of  the  TNO-array.  The  results  are 
also  shown  in  Figure  5.  Although  these  values  tend  to  be  greater  than  the  originally  measured 
ones,  they  still  do  not  match  the  simulation  results.  However,  one  has  to  keep  in  mind  that  only  an 
estimate  could  be  presented,  because  the  hole  diameter  was  measured  from  the  impact  pattern  in 
the  fragment  catcher  in  stead  of  the  real  target  plates.  Moreover,  the  impact  patterns  we  used  only 
showed  the  positions  of  the  craters  with  a  volume  larger  than  0.5  ml  [2],  In  the  simulations  all 
holes  were  taken  into  account,  no  matter  how  small  they  were. 

The  experiments  showed  that  the  hole  diameters  increase  linearly  until  the  sixth  plate  of  the  TNO- 
array  after  which  a  decrease  sets  in.  The  decrease  is  only  indicative  as  it  is  based  on  one  shot.  The 
simulation  shows  an  increase  until  the  seventh  plate  after  which  the  hole  diameters  decrease.  The 
simulation  results  have  also  been  compared  to  the  results  of  experiments  performed  by 
Rheinmetall  and  BWB  against  target  arrays  like  the  TNO-array  [8].  In  general,  the  tendencies 


predicted  by  the  simulations  are  in  agreement  with  the  experimental  results  obtained  by  TNO, 
BWB  and  Rheinmetall.  The  maximum  penetration  capacity  (i.e.  the  number  of  the  last  plate  that 
was  perforated)  can  be  predicted  well.  Both  the  experiments  and  the  simulations  proved  that  up  to 
ten  or  eleven  plates  of  the  target  arrays  can  be  perforated. 


Mass  distribution 

Radiographs  were  used  to  measure  the  fragment  masses  after  a  specific  target  plate.  The  total  mass 
of  break-up  material  and  the  Mott  parameter  were  calculated  from  these  fragment  masses.  The 
relative  mass  of  break-up  material  is  defined  as  the  total  mass  relative  to  the  original  mass  of  the 
projectile's  frangible  part,  and  is  shown  in  Figure  6. 


Plate  no. 


Figure  6  Experimental  and  theoretical  relative  mass  of  break-up  material  within  the  TNO- 
array  (  ★  experiments,  X  simulations,  A  simulations  in  which  only  fragments 
heavier  than  0.01  g  were  taken  into  account) 


The  simulation  results  underestimate  the  experimental  results.  The  experiments  indicated  a  mass 
loss  of  20  percent  at  each  target  plate,  the  simulations  a  mass  loss  of  12  to  15  percent.  The  lost 
mass  represents  the  fragments  that  were  stopped  by  the  target  plates.  However,  in  the  experiments 
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the  lost  mass  also  represents  the  fragments  that  were  too  small  to  measure  (<  0.01  g).  An  adapted 
theoretical  relative  mass  is  shown  in  Figure  6,  which  represents  only  the  fragments  heavier  than 
0.01  g.  Although  the  results  tend  to  be  better  this  way,  there  still  remains  a  gap  between  the 
experimental  and  theoretical  results.  Two  reasons  might  justify  this  gap.  Firstly,  a  limited  number 
of  experiments  was  performed,  and  the  measurement  of  the  fragment  masses  is  not  quite  accurate 
(note  the  relative  mass  of  1 13  %).  Secondly,  the  simulation  results  can  be  tuned  by  decreasing  the 
number  of  particles  into  which  a  fragment  breaks  up  (Nq).  This  way  the  simulation  results  can  be 
improved  with  regard  to  the  relative  mass,  but  the  agreement  with  other  parameters  will  become 
worse.  We  have  chosen  the  number  of  particles  into  which  a  fragment  breaks  up  in  such  way,  that 
the  maximum  penetration  capacity  of  the  fragment  beam  according  to  the  simulation  results 
agrees  with  the  experimentally  found  maximum  penetration  capacity. 

The  mean  fragment  mass  was  determined  from  the  experimental  results  through  the 
determination  of  the  Mott  parameter.  If  the  fragment  mass  distribution  can  be  approximated  by  a 
Mott  distribution,  the  mean  fragment  mass  is  equal  to  two  times  the  Mott  parameter.  Figure  7 
shows  the  Mott  parameter. 


Figure  7 


Experimental  and  theoretical  values  of  the  Mott  parameter  versus  plate  number  of 
the  TNO-array  (  ★  experiments,  X  simulations,  A  simulations  in  which  only 
fragments  heavier  than  0. 1  g  were  taken  into  account) 
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The  same  reasoning  held  at  the  discussion  of  the  relative  mass  holds  for  the  discussion  of  the  Mott 
parameter.  The  deviation  between  the  experimental  and  theoretical  Mott  parameters  depends  on 
the  accuracy  of  the  measurements,  and  on  the  number  of  particles  into  which  a  fragment  breaks 
up.  To  demonstrate  the  influence  of  the  accuracy  of  the  measurements,  another  estimate  of  the 
Mott  parameter  is  shown  in  Figure  7.  Based  on  the  same  fragment  mass  distribution,  we  have 
increased  the  lower  limit  from  0.01  g  to  0.1  g.  The  resulting  Mott  parameter  increases 
dramatically,  especially  after  the  fourth  target  plate  where  lots  of  tiny  fragments  are  present. 

Spatial  distribution 

In  Figure  8,  the  experimental  and  theoretical  cone  angles  are  shown.  The  cone  angle  is  defined  as 
the  aperture  angle  of  the  cone  which  encloses  all  fragments  of  the  fragment  beam. 


1  23456789  10 

Plate  no. 


Figure  8  Experimental  and  theoretical  cone  angles  within  the  TNO-array 
(  it  experiments,  X  simulations) 

Experiments  indicated  that  the  cone  angle  remains  almost  constant  in  the  target  arrays.  Moreover, 
the  cone  angle  seemed  to  be  independent  of  the  target  array  composition.  The  simulation  results 
are  in  good  agreement  with  the  experimental  ones,  although  the  simulations  predict  larger  values. 
The  measured  cone  angles  can  be  underestimated,  since  they  were  measured  from  radiographs. 
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These  radiographs  do  not  show  fragments  smaller  than  0.01  g.  If  the  fragments  smaller  than 
0.01  g  are  left  out,  the  theoretical  cone  angles  will  become  about  10  percent  smaller. 

The  impact  patterns  in  the  fragment  catchers  located  behind  the  target  arrays  were  used  to 
determine  the  number  of  craters,  or  properly  speaking  the  number  of  fragments  and  their  spatial 
distribution.  The  craters  were  divided  into  two  groups;  the  large  craters  with  volumes  larger  than 
0.5  ml,  and  the  small  craters  with  volumes  smaller  than  0.5  ml.  An  equivalent  distinction  has  been 
applied  to  the  simulation  results;  large  fragments  with  a  mass  larger  than  1.786  g,  and  small 
fragments  with  a  mass  smaller  than  1.786  g.  The  relation  between  crater  volume  and  fragment 
mass  was  derived  in  [2].  The  total  number  of  small  and  large  fragments  is  presented  in  Figure  9. 


Figure  9  Experimental  and  theoretical  total  number  of  small  and  large  fragments  within  the 
TNO-array  (small  fragments:  -k  experiments,  X  simulations  1 1  large  fragments: 
▲  experiments,  A  simulations) 

The  number  of  small  and  large  fragments  is  in  the  same  order  of  magnitude  for  the  experimental 
and  simulation  results.  At  the  beginning  of  the  target  array,  the  numbers  agree  well;  moving 
further  in  the  target  plate  array,  the  differences  increase  which  might  be  explained  by  the  large 
amount  of  small  fragments.  The  smaller  the  fragments,  the  worse  the  determination  of  its  number. 
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Kinetic  energy 

The  crater  volumes  were  also  used  to  determine  the  kinetic  energy  of  the  fragment  beam.  It  is 
assumed  that  the  volume  is  proportional  to  the  kinetic  entrgy.  The  total  crater  volume  and  the 
theoretical  relative  total  energy  are  shown  in  Figure  10.  In  contrast  to  the  other  results,  the  crater 
volume  and  the  kinetic  energy  include  the  contribution  of  the  nose  and  rear  part  of  the  FAPDS 
projectile. 


Figure  10  Total  crater  volume  versus  plate  number  as  measured  by  the  experiments  (  *  ), 
and  relative  kinetic  energy  of  the  beam  within  the  TNO-array  as  calculated  by  the 
simulations  (X) 

Although  the  dimensions  of  the  series  in  Figure  10  are  different,  the  tendencies  are  in  perfect 
harmony.  The  decrease  in  energy  is  about  16  to  18  %  per  target  plate.  Therefore,  it  is  obvious  that 
the  maximum  penetration  capacity  is  consistent  for  the  experiments  and  the  simulations. 

An  essential  distinction  between  the  experimental  and  the  simulation  results  exists,  however.  The 
decrease  in  energy  is  equal  to  the  decrease  of  the  measured  fragment  mass  per  target  plate. 
Meaning,  the  velocity  decay  of  the  fragments  is  negligible  according  to  the  experiments.  According 
to  the  simulation  results,  the  decrease  in  energy  is  larger  than  the  decrease  in  fragment  mass,  and 
as  a  result  the  velocity  decay  of  the  fragments  is  not  negligible. 
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4.2  Influence  of  target  composition  and  rifling  twist 

Next  to  the  experiments  with  TNO-arrays,  modified  TNO-arrays  were  used  as  well.  We  will 
discuss  the  results  of  experiments  and  simulations  in  which  one  of  the  aluminium  plates  of  the 
TNO-array  was  replaced  by  a  4  mm  steel  52  plate.  The  steel  plates  were  located  at  position  two 
and  six  respectively.  The  damage  to  the  modified  arrays  is  shown  in  Figure  1 1 . 


Figure  1 1  Experimental  and  theoretical  hole  diameters  of  the  TNO-array  with  steel  plates  at 
position  2  and  6  (steel  plate  at  position  2:  ★  experiments,  X  simulations  I  I  steel 
plate  at  position  6:  ▲  experiments,  A  simulations) 

The  experimental  and  simulation  results  are  in  good  agreement.  Replacement  of  the  aluminium 
plate  by  a  steel  plate  at  position  2  has  hardly  any  influence  on  the  damage  pattern.  Replacement  of 
the  plate  at  position  6  has  a  distinct  influence;  the  hole  diameter  of  the  sixth  plate  decreases 
radically,  and  the  subsequent  hole  diameters  decrease  too.  This  effect  can  be  explained  by  the  so- 
called  overload  concept  [2];  at  the  beginning  of  the  array  the  fragment  beam  has  enough  kinetic 
energy  to  perforate  any  target  plate,  further  down  the  fragment  beam  does  not  have  such 
“overload”  capacity  any  more. 


TNO-report 


PML  23249 1308 


Page 

26 


The  influence  of  the  rifling  twist  on  the  simulation  results  is  clearly  shown  in  Figure  12.  By 
doubling  the  rifling  twist  from  6  to  12  degrees,  a  similar  increase  in  cone  angle  can  be  noticed.  A 
plausible  explanation  has  been  presented  in  the  discussion  of  the  fragmentation  model. 
Experiments  have  demonstrated  the  same  effect  [6]. 


123456789  10 

Plate  no. 


Figure  12  Experimental  and  theoretical  cone  angles  within  the  TNO-array 

(  ★  experiments  with  RT=6",  X  simulations  with  RT=6\  A  simulations  with 
RT=12") 


4.3  Influence  of  parameters  of  the  simulation  model  on  the  results 

The  simulation  model  discussed  so  far  has  been  optimised  to  produce  results,  that  are  consistent 
with  the  experimental  results.  In  this  section,  the  influence  of  mass  distribution  and  scattering 
angle  on  the  simulation  results  is  discussed.  This  way,  we  hope  to  demonstrate  the  urgent  need  for 
more  experimental  data  relevant  to  the  fragmentation  process. 


The  influence  of  the  maximum  number  of  particles  into  which  a  fragment  breaks  up  (Nq)  is 
shown  in  Figure  13  (using  a  Mott  distribution  in  both  situations). 
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Figure  13  Experimental  and  theoretical  hole  diameters  of  the  TNO-array 

(  ★  experiments,  X  Mott  distribution  and  10  <  Nq  <  100,  A  Mott  distribution  and 
10<Nq<20,  O  random  distribution  and  10  <  Nq  <  100;  Nq  is  the  number  of 
particles  into  which  a  fragment  breaks  up) 

As  stated  before,  experiments  indicated  that  the  maximum  number  of  particles  into  which  a 
fragment  breaks  up  should  be  about  20  [2],  The  simulation  model  gives  better  overall  results  when 
the  maximum  number  of  particles  is  increased  to  100.  By  increasing  the  maximum  number  of 
particles,  the  average  mass  of  the  particles  is  being  reduced,  and  as  a  result,  the  penetration 
capacity  of  the  fragment  beam  too.  The  fragment  beam  becomes  wider,  which  results  in  larger 
hole  diameters  as  shown  in  Figure  13. 

The  maximum  number  of  100  (or  even  more)  as  indicated  by  the  simulation  results  is  reasonable. 
The  experimental  values  were  obtained  by  radiographs.  The  particles  on  these  radiographs  are 
hard  to  distinguish.  Moreover,  particles  less  than  0.01  g  are  not  shown  on  these  radiographs. 
Therefore,  the  experimental  values  are  underestimated. 

The  influence  of  the  type  of  mass  distribution  is  shown  in  Figure  1 3  as  well.  A  Mott  distribution  is 
compared  to  a  random  distribution  (using  100  as  the  maximum  number  of  particles  into  which  a 
fragment  breaks  up  in  both  situations).  The  random  distribution  consists  of  a  few  larger  particles 
and  lots  of  tiny  particles.  As  a  result,  even  after  the  perforation  of  several  plates,  a  few  heavy 
particles  will  remain.  The  penetration  capacity  is,  amongst  other  things,  related  to  the  existence  of 
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these  heavy  particles.  A  Mott  distribution  gives  a  somewhat  more  uniform  mass  distribution  with 
regard  to  the  random  distribution.  Therefore,  the  penetration  capacity  will  be  less  and  the 
fragment  beam  wider  whenever  a  Mott  distribution  instead  of  a  random  distribution  is  used. 

As  mentioned  above,  an  increase  in  the  number  of  particles  into  which  a  fragment  breaks  up  (Nq) 
results  in  a  wider  fragment  beam.  Another  possibility  to  control  the  fragment  beam  width  is  to 
change  the  scattering  angle  for  the  individual  fragments  (04)  as  shown  in  Figure  14. 


Figure  14  Experimental  and  theoretical  cone  angles  within  the  TNO-array 

(  *  experiments,  X  Mott  distribution  and  10  <  No  <  100  and  ap2*,  A  Mott 
distribution  and  10<No<100  and  04=3%  O  Mott  distribution  and  10<Nq<20  and 
ctj=3") 

The  scattering  angle  for  the  individual  fragments  does  not  only  affect  the  fragment  beam  width, 
but  also  the  penetration  capacity  of  the  beam.  The  scattering  angle  of  2  degrees  in  combination 
with  the  maximum  number  of  particles  of  100  proved  to  give  the  best  overall  results.  In  Figure  14 
the  simulation  results  using  the  original  experimental  values  are  shown  too.  For  this  particular 
situation  the  fragment  beam  width  is  not  that  bad,  but  the  maximum  penetration  capacity  is  far 
out  of  line. 
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Concluding  remarks 

During  the  development  of  the  simulation  model,  we  tried  to  match  the  simulation  results  to  the 
experimental  ones.  Since  dependencies  are  incorporated,  this  matching  process  is  quite 
complicated.  For  example,  if  we  decrease  the  scattering  angle  for  the  individual  fragments,  the 
cone  angle  decreases,  but  at  the  same  time  the  penetration  capacity  of  the  fragment  beam 
increases.  So,  on  the  one  hand  we  achieve  better  results  for  the  cone  angle,  but  on  the  other  hand 
worse  results  for  the  hole  diameters.  From  all  parameters  which  could  be  compared,  the 
penetration  capacity  of  the  fragment  beam  has  been  marked  as  the  most  important  parameter. 

The  simulation  model  probably  can  be  improved,  if  more  information  on  the  fragmentation 
process  of  the  individual  fragments  is  available.  Planned  experiments  with  frangible  FSPs 
hopefully  will  provide  this  information. 

5  THE  IMPLEMENTATION  OF  THE  FAPDS  SIMULATION  MODEL  IN  A 

VULNERABILITY  ASSESSMENT  CODE 

Because  of  the  encouraging  simulation  results,  it  was  decided  to  implement  the  FAPDS  simulation 
model  in  a  vulnerability  assessment  code.  The  simulation  model  proves  that  the  break-up 
mechanism  of  the  individual  fragments  and  the  effect  of  parasitism  can  be  implemented  in  a 
straightforward  way  in  a  shot  line  based  vulnerability  code.  However,  the  final  implementation 
gives  rise  to  some  as  yet  unresolved  problems. 

The  major  problem  is  to  find  a  way  to  deal  with  synergistic  effects.  In  the  TNO-array,  large  holes 
were  created,  and  cracks  were  formed.  The  target  plates  were  often  bent  and  petals  were  formed. 
Little  is  known  about  these  phenomena.  Models  to  describe  these  mechanisms  in  order  to 
implement  them  in  a  terminal  ballistics  model  still  have  to  be  developed. 

Another  major  problem  is  the  difference  between  real  targets  and  the  simple  target  plate  arrays. 
Live  firings  proved  that  for  instance  electronic  compartments  and  fuel  tanks  have  a  significant 
influence  on  the  damage  pattern.  Replacing  the  electronic  compartments  and  fuel  tanks  with 
simple  boxes  is  not  sufficient  to  obtain  reliable  simulation  results. 

The  FAPDS  simulation  model  is  based  on  a  shot  line  method;  for  every  fragment/particle  a  shot 
line  is  created.  Along  these  shot  lines  calculations  are  performed  to  determine  the  extent  of 
penetration.  Using  a  simple  target  array  with  10  plates  will  result  in  about  100,000  shot  lines  to  be 
evaluated.  Because  of  the  Monte  Carlo  character  of  the  simulation  model,  several  simulation  runs 
(about  50)  have  to  be  performed  to  obtain  statistically  significant  results. 
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One  possible  solution  to  reduce  the  number  of  shot  lines  to  be  evaluated  is  to  define  thresholds. 
For  instance,  the  minimum  mass  necessary  to  perforate  the  weakest  part  of  the  target  could  be 
calculated.  Then  the  shot  lines  for  fragments  with  masses  less  than  the  minimum  mass  do  not  have 
to  be  evaluated  in  a  simulation  run.  This  method  has  already  been  applied  to  the  results  as 
presented  in  Chapter  4.  We  have  defined  a  threshold  of  0.001  g,  which  represents  the  minimum 
mass  necessary  to  perforate  the  2  mm  aluminium  plates  of  the  TNO-array. 

Another  approach  which  could  be  considered  is  the  fragment  density  method.  With  this  method, 
individual  fragments  are  no  longer  traced  through  the  target,  but  they  are  combined  into  one 
single  beam.  The  beam  has  to  be  characterized  by  at  least  two  parameters,  the  beam  width  and  the 
penetration  capacity  of  the  beam.  Effects  like  breaking  up,  parasitism,  and  synergism  do  not  have 
to  be  described  in  detail  for  the  individual  fragments.  Using  this  approach,  one  has  to  tackle  a  new 
problem,  how  to  determine  the  characteristic  parameters  of  the  beam.  A  preliminary  study  [9] 
indicated  that  the  influence  of  the  beam  width  on  effectiveness  computations  cannot  be  neglected. 
The  influence  depends  on  the  component  dimensions  of  the  target  model  used  in  the  effectiveness 
calculations. 

Both  methods  (shot  line  and  beam)  have  advantages  and  disadvantages.  At  the  moment,  we  are 
examining  which  method  can  be  applied  best  for  vulnerability  assessment  studies.  Probably  a 
mixture  of  both  methods,  in  which  a  limited  number  of  shot  lines  will  be  evaluated,  will  be  chosen. 
The  shot  lines  represent  the  fragment  beam.  Penetration  equations  will  be  applied  for  every  shot 
line  (just  like  the  THOR  equations).  Afterwards,  the  penetration  data  of  the  individual  shot  lines 
will  be  combined  to  obtain  the  final  damage  inflicted  upon  the  target. 
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6  CONCLUSIONS 

A  semi-empirical  terminal  ballistics  model  for  FAPDS  munition  has  been  developed.  This  model 
is  based  on  relatively  simple  descriptions  for  the  penetration  and  fragmentation  processes.  The 
effect  of  parasitism  is  included.  The  effect  of  synergism  has  been  included  in  the  aggregation  of 
results. 

Comparison  of  simulation  results  with  experiments  for  35  mm  FAPDS  munition  against  target 
plate  arrays  is  satisfactory.  Two  major  characteristics,  i.e.  the  maximum  penetration  capacity  and 
the  fragment  beam  width,  can  be  calculated  accurately  with  aid  of  the  simulation  program.  The 
influence  of  the  target  array  composition  and  the  rifling  twist  of  the  gun  on  the  simulation  results 
can  be  predicted  well. 

Further  improvement  of  the  simulation  results  is  possible,  when  more  information  on  the  details  of 
the  fragmentation  process  becomes  available.  It  is  expected  that  experiments  with  frangible  FSPs 
will  provide  this  information. 

One  has  to  keep  in  mind  that  the  results  presented  are  valid  for  the  35  mm  FAPDS  projectile  as 
described  in  Chapter  2  only.  Other  kinds  of  FAPDS  projectiles  might  require  new  or  adapted 
models  to  describe  the  fragmentation  process. 

Although  the  problems  affiliated  with  the  implementation  of  the  terminal  ballistics  model  in  a 
vulnerability  code  have  not  been  resolved  yet,  it  is  expected  that  these  problems  will  be  solved  in 
the  near  future. 
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9  LIST  OF  SYMBOLS 


A 

average  presented  area 

m“ 

C 

projectile  calibre 

mm 

D 

projectile  diameter 

mm 

L 

line-of-sight  thickness 

mm 

m 

fragment  mass 

g 

m 

mean  fragment  mass 

g 

min 

fragment  mass  before  penetration  of  a  target 

g 

mout 

N(m) 

N0 

residual  fragment  mass  after  penetration  of  a  target 

cumulative  number  of  fragments  with  mass  >  m 

maximum  number  of  particles  into  which  a  fragment  breaks  up 

g 

r 

hole  radius 

mm 

RT 

rifling  twist 

degrees 

S 

shape  factor 

m2  /  kg2/3 

t 

target  thickness 

mm 

V 

fragment  velocity 

m/s 

v0 

muzzle  velocity 

m/s 

vin 

fragment  velocity  before  penetration  of  a  target 

m/s 

vout 

residual  fragment  velocity  after  penetration  of  a  target 

m/s 

vt  max 

maximum  tangential  fragment  velocity 

m/s 

a 

scattering  angle 

degrees 

“max 

maximum  scattering  angle 

degrees 

0 

angle  between  the  surface  normal  and  the  fragment  trajectory 

degrees 

<d 

spin  rate 

1/s 

At 

perforation  time 

s 
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ANNEX  A  THE  MODIFIED  THOR  EQUATIONS 


The  THOR  equations  probably  are  the  most  familiar  equations  to  calculate  the  residual  mass  and 
velocity  of  fragments  after  penetration  of  a  target  plate.  The  THOR  equations  were  derived 
experimentally  from  ballistic  tests  [7}  for  a  large  number  of  combinations  of  compact  steel 
fragments,  target  materials  and  set-ups. 

The  THOR  equations  are  valid  for  homogeneous  "infinite"  target  plates.  As  explained  in 
Chapter  3,  this  specific  condition  is  not  always  valid  for  FAPDS  munition.  Situations  occur  in 
which  fragments  encounter  holes  in  the  target  plates.  Depending  on  the  dimensions  of  the 
fragments  and  of  the  holes,  some  fragments  can  fly  through  the  holes  without  any  interaction  with 
the  target  plate  (parasitism).  Other  fragments,  that  are  larger  than  the  holes,  will  interact  with  the 
target  plate,  and  as  a  result  will  break  up.  Besides,  the  energy  of  the  fragments  will  be  reduced,  and 
the  holes  in  the  target  plate  will  be  enlarged.  To  calculate  the  residual  mass  and  velocity  of  these 
semi-parasitic  fragments,  we  have  derived  the  so-called  modified  THOR  equations.  The  derivation 
of  these  equations  is  discussed  below. 


The  THOR  equations  are  the  basics  for  the  modified  THOR  equations.  Therefore,  we  will 
summarize  some  part  of  Chapter  3  in  which  the  THOR  equations  were  discussed. 

Consider  the  configuration  shown  in  Figure  A.  1 . 


mout,vout 


Figure  A.  1  A  fragment  penetrating  a  target 

A  fragment  with  mass  mjn  and  velocity  Vjn  hits  the  target.  The  angle  between  the  velocity  vector  of 
the  fragment  and  the  normal  of  the  target  surface  is  defined  as  9.  During  the  penetration,  the 
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velocity  of  the  fragment  as  well  as  its  mass  reduces.  If  the  fragment  strikes  the  target  with  sufficient 
energy,  it  will  perforate  the  target.  The  residual  mass  mout  and  velocity  vout  can  be  calculated 
with  aid  of  the  following  THOR  equations: 


10CI  *  (t*A)C2  *  minC3  * 


C4  *  C5 


mout  =  min  -  1 0K 1  *  (t*A)K2  *  minK3  *  K4  *  vinK5  (A.  2) 

where  C 1  ..C5  and  K1  ..K5  are  material  constants 

t  =  target  thickness 

A  is  the  average  presented  area  of  the  fragment  calculated  from  the  shape  factor  S 
and  the  mass  mjn: 

A  =  S  *  min2/3  (A.  3) 

If  the  residual  mass  or  velocity  is  less  than  or  equal  to  zero,  the  fragment  is  not  capable  of 
perforating  the  target. 

For  the  derivation  of  the  modified  THOR  equations,  we  use  the  situation  as  shown  in  Figure  A. 2. 


Figure  A.2  A  fragment  flying  through  a  hole  in  the  target 
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A  fragment  with  mass  mjn  and  velocity  Vjn  hits  the  target  at  the  location  where  a  hole  is  present. 
The  area  of  the  hole  is  A^0je  (frontal  area).  The  presented  area  of  the  fragment  is  defined  by 
equation  A.  3. 

We  will  only  consider  the  situations  of  normal  impact,  which  are  of  most  interest  for  FAPDS 
munition  (the  striking  angles  are  small,  up  to  about  3‘).  As  a  result,  the  modified  THOR  equations 
will  be  valid  for  normal  impacts  and  for  small  impact  angles.  We  assume  that  the  centre  of  mass  of 
the  fragment  flies  through  the  centre  of  the  hole.  This  way,  we  neglect  troublesome  phenomena, 
such  as  tumbling  of  the  fragment  which  might  occur  if  the  fragment  strikes  the  edge  of  the  hole. 
Finally,  we  assume  that  the  cross  section  of  the  hole  is  equal  along  the  axis  of  the  hole. 

Regarding  the  situation  shown  in  Figure  A.2,  two  possibilities  arise: 

1  The  fragment  is  smaller  than  the  hole. 


A  <  Ajiole 


(A.4) 


or 

2  The  fragment  is  larger  than  the  hole,  or  has  an  equal  dimension. 


A  >=  Ahole  (A. 5) 

The  first  possibility  is  quite  simple  to  deal  with;  the  fragment  flies  through  the  hole  without  loss  of 
energy.  Meaning: 


mout  hole  ~  min  >  anc*  vout  hole  -  vin  (A- 6) 

where  mout  hole  and  vout  hole  represent  the  residual  mass  and  velocity  of  a  fragment  that  went 
though  a  hole  in  the  target  plate 

The  second  possibility  is  far  more  complex.  The  energy  of  the  fragment  will  be  reduced,  because 
the  fragment  has  to  penetrate  the  target.  The  residual  mass  and  velocity  of  the  fragment  will  be 
larger  with  regard  to  the  situation  in  which  no  hole  is  present  in  the  target  plate.  As  a  result,  the 
existing  hole  will  be  enlarged. 
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Consider  the  situation  in  which  no  hoie  is  present  in  the  target  plate.  The  energy  of  the  fragment 
before  and  after  penetration  is  defined  by  equations  A. 7  and  A. 8. 


m;, 


Mn 


(A. 7) 


^out  ~  2 


m 


out  vout 


(A. 8) 


The  energy  Eperf,  i.e.  the  energy  needed  to  perforate  the  target  plate,  is  defined  as  the  difference. 


Eperf  -  ^in  '  ^out  (A. 9) 

The  derivation  of  the  modified  THOR  equations  is  based  on  energy  conservation  laws.  What  will 
be  done  is  to  add  the  minimum  amount  of  energy  needed  to  perforate  the  target  plate  (and  create 
a  hole  with  area  Aj,0ie)  to  the  residual  energy  of  the  fragment  if  no  hole  was  present  in  the  target 
plate  (equation  A.8).  From  the  resulting  residual  energy  the  residual  mass  mout  ^ole  and  velocity 
vout  hole  can  *5e  calculated. 

The  minimum  energy  needed  to  create  a  hole  (with  frontal  area  Aj,0je  in  a  target  plate  of  thick¬ 
ness  t)  results  from  equation  A.  10. 


^  Epertfmin>  vin)  _  ^  (  ^in  '  ^out ) 
'd  vin  3  min  d  vin  d  min  " 


(A.  10) 


The  hole  size  is  chosen  to  be  equal  to  the  presented  area  of  the  fragment,  that  created  the  hole.  As 
a  result,  m,n  is  set  by  equation  A.  3 1 . 


min  =  <- 


^hole  , 
S  ' 


3/2 


(A.  11) 
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Equation  A.  10  can  be  solved  by  use  of  equation  A.  1 1  and  the  THOR  equations.  This  is  a  rather 
complicated  problem,  while  the  THOR  equations  are  interrelated,  and  can  only  be  used  if  mout 
and  vout  are  greater  or  equal  to  zero. 

Therefore,  another  approach  has  been  selected  in  which  the  residual  energy  Eout  has  been  set  to 
zero.  Now,  the  THOR  equations  can  be  used  to  calculate  the  value  of  vjn  corresponding  to  the 
minimum  perforation  energy.  This  value  of  vjn  will  be  referred  to  as  vp.  The  value  of  mjn 
corresponding  to  the  minimum  perforation  energy  will  be  referred  to  as  mp.  The  minimum 
perforation  energy  Eperf  min  can  calculated  with  aid  of  equation  A.  12. 


Eperf  min 


(A.  12) 


where 


Ahole 
mp  =  (  5 


3/2 


and  vp  =  minimum  of 


(A.  13) 


10C1  *  (t*Ahole)C2  *  mpC3 


±_\ 

cos  0^ 


C4 


1-C5 


and 


10'K1  *  (t*Aj10je)_K2  *  mpl~K3 


1  ^ 


cos  8 


-K4 


) 


K5 


(A.  14) 


As  said  before,  the  minimum  perforation  energy  (Eperf  min)  *s  to  he  a<*ded t0  the  residual  energy 
of  the  fragment  (Eout)  if  no  hole  was  present  in  the  target  plate.  The  sum  of  those  energies  is  the 
residual  energy  of  the  fragment  (Eout  hole)>  that  penetrated  the  target  plate  with  the  hole. 
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Once  again,  we  have  to  distinct  the  following  three  situations: 

1  The  energy  of  the  fragment  is  sufficient  to  penetrate  even  an  equivalent  target  plate  without  any 
hole. 

Eout  >  0  (A. 1 5) 

2  The  energy  of  the  fragment  is  sufficient  to  penetrate  the  target  plate  while  a  hole  in  the  target 
plate  is  present. 


Eout  <_  0 


and 


Eout  hole  Eout  ■*"  Eperf  min  )  >  ®  (A.  16) 

3  The  energy  of  the  fragment  is  not  sufficient  to  penetrate  the  target  plate,  even  if  a  hole  is 
present  in  the  target  plate. 


Eout  <=  0 


and 


Eout  hole  Eout  +  Eperf  min  )  <_  ®  (A.  17) 

Negative  energies  have  been  used  to  indicate  that  the  fragment  is  not  capable  of  perforating  the 
target  plate. 

For  each  of  the  three  situations  mentioned  above,  the  parameters  mout  hole  and  vout  hole  can  Ee 
derived  form  equation  A.  18. 


Eout  hole  2  *  mout  hole  *  vout  hole^  ~  Eout  +  Eperf 


min 


(A.  18) 
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However,  this  is  only  one  equation,  while  we  have  to  determine  two  parameters.  We  use  the 
THOR  equations  to  solve  this  problem.  As  shown  by  these  equations,  the  residual  mass  and 
velocity  of  a  fragment  is  determined  by  (t*A)^2  and  (t*A)^2  respectively.  The  id.  (t*A)  is  the 
volume  of  target  material,  that  has  to  be  removed  from  the  target  plate  during  penetration.  In 
other  words,  the  id.  (t*A)  represents  some  kind  of  ballistic  resistance,  which  has  to  be  broken. 
Equation  A.  18  can  be  solved  using  this  information.  The  penetration  of  a  target  plate  with  a  hole 
is  compared  to  the  penetration  of  a  target  plate  of  less  ballistic  resistance.  The  ballistic  resistance, 
in  particular  the  id.  (t*A),  will  be  referred  to  as  BRC  (Ballistic  Resistance  Coefficient). 

The  following  equations  can  be  derived  from  the  THOR  equations  at  given  mjn,  v;n  and  0: 


mout  '  (  BRC  )  K2 


(A.  19) 


Vout'(BRC)C2 


(A.  20) 


As  a  result: 


Eout  '  (  BRC  )  K2+2*C2  (A.21) 

Comparing  the  penetration  of  a  target  plate  with  a  hole  to  the  penetration  of  a  target  plate  of  less 
ballistic  resistance  (BRC^ole)’  we  arr*ve  at  die  following  equations: 


mout  hole 


f 


BCfoole 


BRC 


K2 


m 


out 


_  f Eout  hole^ 


K2 


-out 


■J 


K2+2*C2 


m 


out 


(A.  22) 


_  fB^-hole^C^  ^ 
vout  hole  ~  ^  BRC  J 


'out 


c  N  C2 

_  f  Eout  hole^  “  * 

“I  Eout  " JK2+2*C2  vout 


(A.23) 


where  BRC  =  the  ballistic  resistance  coefficient  of  the  target  plate  without  holes. 
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Equations  A. 22  and  A. 23  form  the  modified  THOR  equations.  They  provide  the  solution  for 
situation  no.l,  which  has  been  described  earlier;  the  energy  of  the  fragment  is  sufficient  to 
penetrate  an  equivalent  target  plate  even  without  any  hole. 

The  solution  of  situation  no.  2  (the  energy  of  the  fragment  is  sufficient  to  penetrate  the  target  plate 
while  a  hole  in  the  target  plate  is  present)  is  not  available  yet.  The  modified  THOR  equations  can 
not  be  used,  because  mout  and/or  vout  are  not  valid  (tnout  <=  0  and/or  vQUl  <=  0).  In  the  FAPDS 
simulation  model,  this  problem  has  been  resolved  by  assuming  that  the  hole  in  the  target  plate  will 
be  enlarged,  but  the  residual  energy  of  the  fragment  will  be  negligible. 

The  solution  of  situation  no.  3  (the  energy  of  the  fragment  is  not  sufficient  to  penetrate  the  target 
plate,  even  if  a  hole  is  present  in  the  target  plate)  is  simply  given  by  equation  A.  24 


mout  hole  0  >  anc*  vout  hole  ~  ®  (A. 24) 

The  explanation  is  obvious;  the  fragment  is  not  capable  of  perforating  the  target  plate.  In  the 
FAPDS  simulation  model,  we  assumed  that  the  fragment  will  not  enlarge  the  hole  in  the  target 
plate. 

Summary 

The  following  assumptions  have  been  made  to  derive  the  modified  THOR  equations: 

1  Normal  impact  of  the  fragments  on  the  target  plate.  However,  the  results  may  also  be  used  for 
small  impact  angles. 

2  Phenomena,  such  as  tumbling  of  the  fragment  which  might  occur  if  the  fragment  strikes  the 
edge  of  the  hole,  are  neglected. 

3  The  cross  section  of  the  hole  is  equal  along  the  axis  of  the  hole. 

4  The  hole  size  is  equal  to  the  presented  area  of  the  fragment,  that  created  the  hole. 

5  The  penetration  of  a  target  plate  with  a  hole  can  be  compared  to  the  penetration  of  a  target 
plate  of  less  ballistic  resistance. 

The  following  situations  can  be  distinguished: 

I  The  fragment  is  smaller  than  the  hole,  meaning: 


A  <  fyiole 
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Results: 

mout  hole  _  min»  an<3  vout  hole  -  vin 

II  The  fragment  is  larger  than  the  hole,  and  the  energy  of  the  fragment  is  sufficient  to  penetrate 
an  equivalent  target  plate  even  without  any  hole. 

A  >=  Afooie 


and 


^out  >  0 


Results: 


mout  hole- 


F-out  +  ^perf 


K2 


mm 


-out 


1k2+2*C2*  mout 


vout  hole  - 


^out  +  Eperf 


C2 


mm 


V. 


-out 


^K2+2*C2  ~  vout 


where  Eperf  mm  =  \  *  (— ^  )  3/2  *  vp2 
and  Vp  =  minimum  of 


«*Ahole)C2  * 


.  Ahole 
(  S 


3/2  C3  * 


1  \ 


cos  8j 


C4 


and 


10-K1  *  (t*Ahole)-K2  * 

V 


(  Ahole  )  3/2  (1-K3)  * 


(  1 


cos  0 


-K4 


T  NO- report 


_ _ ANNEX  A 

PML  232491  308  Page 

10 


III  The  fragment  is  larger  than  the  hole  and  the  energy  of  the  fragment  is  sufficient  to  penetrate 
the  target  plate  while  a  hole  in  the  target  plate  is  present. 

A  >=  Ah0ie 


and 


-out 


i 


and 


F'out  hole  (  Eout  +  Eperf  mjn  )  >  0 


Results:  Not  known. 

In  the  FAPDS  simulation  model,  we  have  assumed: 
mout  hole  =  0  >  and  vout  hole  =  0  >  and  the  hole  is  enlarged 

IV  The  fragment  is  larger  than  the  hole,  and  the  energy  of  the  fragment  is  not  sufficient  to 
penetrate  the  target  plate,  even  if  a  hole  is  present  in  the  target  plate. 


A >  Ah0ie 


and 


^out  <_  0 


and 


Results: 


mout  hole  ~  0  >  and  vout  hole  -  ® 


TNO-report 


_  ANNEX B 

PML  232491308  Page 

1 


ANNEX  B  THE  FAPDS  SIMULATION  COMPUTER  PROGRAM 

The  FAPDS  simulation  computer  program  has  been  described  in  Chapter  3.  The  structure  of  the 
program  is  discussed  in  this  Annex. 

The  FAPDS  simulation  computer  program  is  built  modularly,  as  illustrated  in  Figure  B.  1 . 


Figure  B.  1  Structure  of  FAPDS  simulation  computer  program. 

The  main  module  entided  'FAST'  uses  the  modules  ’/m'r'and  'Simulate'.  'Simulate'  contains  the 
simulation  loop  in  which  the  'Evaluate  Event'  module,  the  key  to  the  simulation,  is  used.  Here  the 
physical  processes  are  simulated  which  control  the  penetrauon  and  breaking  up  of  individual 
fragments.  To  that  end  the  module  uses  a  number  of  other  modules.  A  short  descripuon  of  the 
modules  is  given  below. 
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FAST: 

The  module  which  contains  the  main  FAPDS  simulation  loop. 

Inti: 

The  module  containing  the  procedures  to  initialize  the  database  containing  the  model  of  the  target 
and  to  visualize  that  model. 

Simulate: 

The  module  containing  the  procedures  to  perform  the  FAPDS  simulation  loop. 

Evaluate  Event: 

This  module  consists  of  two  modules:  'Evaluate'  and  'Event  list'.  The  latter  contains  the  procedures 
to  maintain  the  event  list.  'Evaluate'  contains  the  procedures  to  evaluate  a  single  event. 

Penetrate: 

The  module  containing  the  procedures  to  calculate  the  parameters  involved  in  the  penetration 
process  whenever  a  fragment  hits  the  target. 

Thor: 

The  module  containing  the  (modified)  THOR  equations  to  calculate  the  residual  fragment  mass 
and  velocity  after  penetration. 

Timelapse: 

The  module  containing  the  procedures  to  calculate  the  time  it  takes  a  fragment  to  perforate  the 
target. 

Breakup: 

The  module  containing  the  procedures  to  break  up  a  fragment  into  pieces,  with  regard  to  the 
masses. 

Scatter: 

The  module  containing  the  procedures  to  scatter  the  pieces  into  which  a  fragment  breaks  up 
spatially. 

Holes: 

The  module  containing  the  procedures  to  deal  with  the  holes  created  in  the  target  by  the 
fragments. 

Shoot  ray: 

The  module  containing  the  procedures  to  calculate  the  intersections  of  the  fragment  trajectories 
with  the  target. 

Most  of  these  modules  are  straight  forward  implementations  of  the  models  discussed  in  Chapter  3. 
The  'Evaluate'  module  is  a  more  complicated  one,  because  each  event  in  the  event  list  has  to  be 
evaluated. 
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The  evaluation  of  a  single  event  involves  the  following  aspects,  as  illustrated  in  Figure  B.2.  The 
first  step  of  the  evaluation  process  is  to  determine  whether  the  fragment  flies  through  a  hole.  If  so, 
and  the  fragment  is  smaller  than  the  hole,  the  mass  and  velocity  of  the  fragment  remain  unaltered 
and  the  fragment  trajectory  (shoot  ray)  is  evaluated  further  until  it  hits  another  part  of  the  target. 
The  collision  time  is  calculated  and  a  new  event  is  added  to  the  event  list.  If  the  fragment  is  larger 
than  the  hole,  the  residual  mass  and  velocity  of  the  fragment  is  calculated  by  use  of  the  modified 
THOR  equations.  If  the  fragment  is  actually  capable  of  perforating  the  target,  the  hole  is  adapted 
and  new  vectors  for  the  resulting  particles  are  created. 

If  no  hole  is  present  at  the  impact  point,  the  next  step  is  to  calculate  whether  the  fragment  is  able 
to  perforate  the  target.  If  the  fragment  is  able  to  perforate  the  target,  the  residual  mass  and  velocity 
are  calculated,  as  well  as  the  location  of  the  exit  point  and  the  required  perforation  time.  At  the 
exit  point  the  fragment  breaks  up  and  new  vectors  (in  fact  trajectories  of  pieces)  are  created.  For 
each  vector  a  possible  new  hit  point  is  determined,  and  a  new  event  is  added  to  the  event  list. 
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